Generation of Bright Two-Color Continuous Variable Entanglement 
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We present the first measurement of squeezed-state entanglement between the twin beams pro- 
duced in an Optical Parametric Oscillator (OPO) operating above threshold. Besides the usual 
squeezing in the intensity difference between the twin beams, we have measured squeezing in the 
sum of phase quadratures. Our scheme enables us to measure such phase anti-correlations between 
fields of different frequencies. In the present measurements, wavelengths differ by » I nm. Entan- 
glement is demonstrated according to the Duan el al. criterion [Phys. Rev. Lett. 84, 2722 (2000)] 
A 2 p_ + A 2 ij+ = 1.47(2) < 2. This experiment opens the way for new potential applications such as 
the transfer of quantum information between different parts of the electromagnetic spectrum. 

PACS numbers: 03.67.-a, 03.67.Mn, 42.50.Dv, 42.50. Yj 



The field of quantum information has recently at- 
tracted great interest, owing to potential applications in 
information storage, communication, and computing 0. 
Entanglement is viewed as a key resource for these ap- 
plications, especially for quantum communication. A va- 
riety of physical systems presenting entanglement have 
been investigated both theoretically and experimentally. 
The vast majority of experiments concentrate on discrete 
variable systems, such as trapped ions J2J, few photon 
electromagnetic fields in cavity QED [3j, spontaneous 
parametric downconversion , and nuclear magnetic res- 
onance 0. On the other hand, in recent experiments, 
continuous variable systems have been studied, includ- 
ing light beams [|| and atomic samples Q. The number 
of such experiments, however, is still relatively small in 
comparison with discrete variable systems. 

The first experimental demonstration of continuous 
variable entanglement used a continuous-wave (CW) 
optical parametric oscillator (OPO) operating below 
threshold The two squeezed vacuum outputs were 
shown to possess Einstein Podolsky Rosen (EPR) type 
correlations. Most recent experiments use a nonlinear 
interaction to produce squeezed fields, which are then 
combined in a beamsplitter to generate entanglement || . 
Conversely, a beamsplitter transformation can also be 
used to generate a squeezed beam from an entangled in- 
put In these experiments, it is mandatory to have 
fields of the same frequency. 

In this Letter, we present the first (to our knowledge) 
measurement of continuous variable entanglement be- 
tween bright fields of truly different frequencies. Even 
before the first experiment 0, it was predicted that 
the above-threshold OPO should produce entangled twin 
beams 0,0. So far, this prediction had not yet been 
verified, owing to the difficulty of measuring phase prop- 
erties of the twin beams. We have succeeded in mea- 
suring quantum correlations between the sum of phase 
quadratures of non-degenerate twin beams. 

Bipartite continuous variable entanglement 0] can be 
tested according to a criterion established by Duan et 



al. and also by Simon 0] . The criterion is based on 
the total variance of EPR- type operators. For operators 
Xi and f>i that obey position-momentum commutation re- 
lations, they consider the variances of combinations such 
as x\ + £2 and p\ — p2- The quadratures of electromag- 
netic fields satisfy such commutation relations. We fo- 
cus here on the so-called amplitude p and phase q field 
quadratures, defined by: 
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where a is the field annihilation operator, (f> is an ar- 
bitrary phase and, for a macroscopic field with a well- 
defined mean amplitude, (q) = 0. In terms of these oper- 
ators, inseparability (entanglement) is demonstrated by 
a violation of the inequality: 
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where the Standard Quantum Level (SQL) is normalized 
to 1 for each combination of quadratures. In order to 
simplify notation, we will refer to {p\ ijb) / as p± , and 
to (qi ± q2)/V2 as q± [thus, A 2 p_ < 1 signals squeezing 
in the difference of amplitude quadratures]. If both p_ 
and g+ are squeezed, inequality (J2J) is violated and we 
have squeezed-state entanglement. 

It is easy to understand why entanglement is expected 
in the CW OPO operating above threshold. The OPO 
consists of a nonlinear crystal inside a resonant cav- 
ity, in which parametric downconversion takes place. The 
cavity feeds back the downconverted fields, leading to 
stimulated parametric gain and hence to an oscillation 
threshold. Since the primary downconversion process in- 
volves creating pairs of photons by annihilation of pump 
photons, one naturally expects strong correlations be- 
tween the intensities of the twin beams: positive in- 
tensity fluctuations of one beam correspond to positive 
intensity fluctuations of the other beam. These corre- 
lations, however, are frequency-dependent. For power 
spectrum analysis frequencies larger than the OPO cav- 
ity bandwidth, correlations tend to disappear since, for 
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times shorter than the cavity lifetime, a photon can exit 
the cavity while its "twin" still remains inside. Squeez- 
ing in the intensity difference was already observed back 
in 1987 ^5|. On the other hand, energy conservation 
(luq = u>i + u>2, where indices 0, 1, and 2 refer to pump, 
signal and idler beams, respectively) and phase match- 
ing imply strong anti-correlations between phase fluctua- 
tions: positive phase fluctuations of one beam correspond 
to negative phase fluctuations of its "twin" . This is ex- 
actly the situation discussed following J5J): entanglement 
occurs if these (anti)correlations lead to fluctuations be- 
low the SQL. However, squeezing in the phase-sum had 
not been measured to date. 

The first prediction of entanglement in the above- 
threshold OPO was made by Reid and Drummond [Io| . 
who even suggested a way of measuring phase fluctu- 
ations and analyzed the low frequency effects of phase 
diffusion (llj . A recent detailed prediction, taking into 
account the effects of pump noise and cavity detunings 
for the three fields, for a triply resonant OPO, was pre- 
sented in [16j. The basic difficulty in measuring (phase) 
quadrature fluctuations, in contrast with intensity fluc- 
tuations, is that the standard technique, homodyne de- 
tection, calls for a local oscillator having a well-defined 
phase relationship with the field to be measured. In the 
OPO, this is difficult to implement, since the frequencies 
of the twin beams are usually different, depending on the 
oscillating modes, and vary from one realization to the 
next. Hence, two local oscillators would be required: one 
for each beam. It would also be necessary to phase-lock 
these fields to the twin beams. 

One way to overcome this difficulty is to force the OPO 
to oscillate in a strictly frequency-degenerate situation. 
This is technically challenging and has been done by two 
groups, using different approaches [ItI Hl| . Our strategy 
is to perform self-homodyne measurements, without the 
use of local oscillat ors, by a frequency-dependent reflec- 
tion of each beam |l9j. If one considers the field as a 
mean value at a carrier frequency with noise sidebands 
at some analysis frequency, a frequency-dependent re- 
flectivity entails different phase shifts for the carrier and 
sidebands. Consequently, different quadrature fluctua- 
tions can be projected onto amplitude fluctuations (with 
respect to the phase of the mean field). A detailed de- 
scription for a single field reflected off an optical cavity 
was given by Galatola et al. |2(1 | . For an imperfect cavity, 
in which vacuum leaks from the outside through the mir- 
rors, the reflected beam amplitude noise power spectrum 
Sn(fl) can be written as 



mission coefficients through the relations: 
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Amplitude reflection r(A) and transmission t(A) coeffi 
cients can be simply written as: 

r 2 exp(iA/du ac ) 
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where A is the detuning between the incident field cen- 
tral frequency and the cavity resonance frequency, 5v ac 
is the cavity bandwidth (FWHM), and Q, is the analysis 
frequency. Cavity input mirror amplitude reflection and 
transmission coefficients are denoted by r\ and t\ , while 
t 2 is defined so that all internal losses A obey the relation 
t\ = A = \-r\. 

For analysis frequencies larger than \[2 8v^ Cl it is pos- 
sible to completely convert incident phase fluctuations 
into amplitude fluctuations of the reflected beam. As 
a matter of fact, if this condition is satisfied, then, for 
A = 0, |(? p | 2 s» 1 and we recover amplitude fluctuations. 
Amplitude fluctuations are also recovered for A 3> 5v ac , 
for which \g p \ 2 — > 1. For A = ±£2/ ac /2, then g q « 1 
and incident phase fluctuations are projected onto am- 
plitude fluctuations of the reflected beam. We use one 
analysis cavity for each beam and scan their frequencies 
synchronously. In this way, we are always measuring the 
same quadrature for each field, with respect to its mean 
value, regardless of the frequency difference between the 
fields. 

Analysis Cavity 1 
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Analysis Cavity 2 



S R {n) = \g p \ 2 S p (Sl) + \g q \ 2 S q (Q) + \g vp \ 2 + \g vq \ 2 , (3) 

where S P (Q) and S q (Q) are the incident beam amplitude 
and phase noise, respectively, and g p , g q , g vp and g vq are 
coefficients that depend on cavity reflection and trans- 



FIG. 1: Sketch of the experimental setup (details explained 
in the text). 

Our experiment is performed with a triply resonant 
CW OPO operating above threshold [2]]]. The pump 
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laser is an ultrastable diode-pumped frequency-doubled 
Nd:YAG source (Innolight Diabolo) at 532 nm, with a 
second output beam at 1064 nm, which we use for align- 
ment purposes. The nonlinear crystal is a 12 mm-long 
type-II High Gray Tracking Resistant KTP from Raicol. 
Crystal temperature is kept near 24°C, with a stability 
of the order of 10 mK, by means of a peltier element. 
The cavity is a quasi-confocal linear Fabry Perot cavity, 
with input mirror reflectivities equal to 89% at 532 nm 
and greater than 99.8% at 1064 nm. Output mirror re- 
flectivities are greater than 99.8% at 532 nm and 95% 
at 1064 nm. The typical threshold power is 60 mW and 
OPO cavity bandwidth is 5v OFO — 52(1) MHz. Noise 
in the difference of signal and idler intensities has a sta- 
ble value A 2 p_ = 0.63(1) registered by the photodetec- 
tors, or —2.0(7) dB. According to [lj|, A 2 g + is signifi- 
cantly affected by pump excess phase noise. Our laser 
presented excess noise up to 25 MHz, which we had to 
filter by transmission through a ring cavity, with a band- 
with 5i/f = 2.3(1) MHz. In this way, the pump beam was 
shot-noise-limited for frequencies above 15 MHz. 

The experimental setup is sketched in fig-^ The pump 
beam at 532 nm is sent through the filter cavity and then 
mode-matched to the OPO cavity. The orthogonally po- 
larized twin beams produced, with powers of the order 
of a few milliwatts each, are separated by a polarizing 
beamsplitter cube (PBS) and each directed to a tunable 
ring analysis cavity. For our working crystal tempera- 
ture, wavelengths of signal and idler beams can differ 
by 0.8 nm to 0.9 nm. Analysis cavity bandwidths are 
dv ac = 14(1) MHz. The reflected field is detected by 
a high quantum efficiency [85(3)%] photodiode (Epitaxx 
ETX 100). The photocurrent is preamplified and the DC 
and High Frequency (HF) components separated. The 
HF components are sent to a demodulating chain, where 
they are mixed with a sinusoidal reference at the analy- 
sis frequency n = 27 MHz (RBW of 600 kHz). As the 
analysis cavities' resonances are swept over time, vari- 
ances of each individual noise component, of their sum, 
and of their difference are calculated. The number N 
of points used, which is proportional to the acquisition 
time, is large enough to guarantee a well defined variance 
and small enough to correspond to a very small change 
in cavities' detuning (this is equivalent to a spectrum an- 
alyzer VBW, in our case approximately 1 kHz). 

Sum and difference noise spectra recorded as functions 
of the synchronous analysis cavity frequency scans are 
presented in fig. A scan over ±4 5v ac is presented in 
fig. El (a). From the sum and difference spectra, one can 
easily notice that the quadrature being measured on each 
beam alternates between amplitude and phase. In par- 
ticular, we recognize the amplitude quadrature behavior 
for detunings |A| > 3 5v ac and for A = 0. One observes 
that A 2 p_ at A = does not recover the squeezing ob- 
served for |A| > 3 Sv ac , owing to a lack of experimental 
resolution (since the cavity frequency scan is large com- 
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FIG. 2: Measurements of sum (red) and difference (blue) of 
quadrature fluctuations as functions of analysis cavities' de- 
tunings. In (a) the detuning spans the range ±4 Sfac- In 
(b), a local scan is presented, in which we clearly observe 
A 2 p_ < 1 for A = and A 2 q+ < 1 for |A| = 0.5 8u ac , char- 
acterizing squeezed-state entanglement. The solid curves are 
theoretical fits to the data (eq|3Jl. 



pared to Si> ac , very few data points are obtained at any 
particular detuning). Phase quadrature is measured for 
|A| = 0.5 8v ac and, for the particular choice of analysis 
frequency we use, also at |A| = 1.65 8v ac . For all other 
detunings we measure a linear combination of amplitude 
and phase quadratures. In fig. [2(b), we present a local 
scan, with increased resolution, over ±1 5v ac . Phase sum 
squeezing is observed, with A 2 q+ = 0.84(2). 

The solid curves in fig. (a) and (b) represent the- 
oretical fits to the data, taking into account the 
transfer of quadrature noise to amplitude noise of the 
fields reflected from the analysis cavities (eq. 0) • Apart 
from scale factors and curve central position, the relevant 
free parameters are the combined quadrature fluctuations 
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A 2 p + , A 2 p_, A 2 g + , and A 2 g_. All other variables re- 
quired in eq. [21 such as cavity bandwidth and analysis 
frequency, are independently measured and employed as 
constants in the fitting, which is done either for the sum 
or for the difference in each part of fig. [21 with equally 
good results. 

In the first spectra we recorded, with pump power 
equivalent to 1.5 times the threshold power, we observed 
str ong excess noise in the phase sum, as was also found 
in [l£| ■ This excess noise is not predicted by the stan- 
dard OPO linearized theory ^(J and it is still not clear 
whether a full quantum theory can account for it. 
Although the calculations in [la] are carried out with the 
linearized theory, they provided us with useful indica- 
tions: squeezing in the phase sum improves for nonzero 
detunings of signal and idler beams with respect to the 
OPO cavity. This behavior is also observed as pump 
power decreases, approaching the threshold power. Con- 
sistently with these indications, we could only observe 
phase sum squeezing very close to threshold (approxi- 
mately seven percent above threshold). In this situation, 
the OPO is very unstable, hindering the measurements. 

In terms of the entanglement criterion of eq. [2J using 
A 2 q+ = 0.84(2), we obtain, A 2 p_ + A 2 q + = 1.47(2) < 2, 
a clear violation of the inequality, characterizing, for the 
first time, entanglement between bright beams of truly dif- 
ferent frequencies. Eq. |2 is a necessary and sufficient 
separability condition for gaussian states, which are pre- 
dicted for the OPO and are consistent with our data. We 
have substantial losses from the output of the OPO to the 
detectors. If the overall detection efficiency is r), the mea- 
sured variances are related to the "true" variances (e.g., 
A 2 p') by A 2 p = rj{A 2 p' — 1) + 1. Correcting for losses, 
equal to 28(2)%, we obtain A 2 p'_+A 2 q' + = 1.26(4). An- 
other, more stringent criterion, is the so-called EPR crite- 
rion, which enables one to infer variances on one beam, as 
functions of variances on the other beam 101. This crite- 



rion implies a violation of the inequality A 2 pj n f 
1, where 
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(pi). An analogous relation holds for 



A <7i n f. Our measured data does not violate the inequal- 
ity, resulting in a product equal to 1.09(4). However, cor- 
recting for losses, the product is A 2 p[ nf A 2 q[ ni = 0.91(5), 
violating this second criterion as well. 

The entanglement generated in this system can be sub- 
stantially improved. The intensity difference squeezing 
measured in [l£| reached the impressive value of 9.7 dB, 
which would lead to even stronger violations of the above 
criteria. One other problem in our system is the unstable 
operation very close to threshold. However, very stable 
operation only a few percent above threshold can be ob- 
tained, as shown in 



In summary, we have demonstrated, for the first time, 
bright two-color squeezed-state continuous variable en- 
tanglement. Applications to quantum information, such 
as quantum key distribution [2j| and quantum telepor- 
tation |24| with continuous variables, can be easily envis- 
aged. Our measurement scheme has interesting proper- 
ties for both. Quantum key distribution with squeezed- 
state entanglement usually requires sending the local os- 
cillator, in addition to the entangled field, so as to en- 
able homodyne measurements, a requirement which is 
not necessary in our case. Teleportation, on the other 
hand, has been restricted to fields of the same frequency. 
Two-color entanglement opens the way for distributing 
quantum information between different parts of the elec- 
tromagnetic spectrum. 
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